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Abstract 

The leading electroweak MSSM two-loop corrections to the electroweak precision 
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observables are calculated. They are obtained by evaluating the two-loop 0(a 
0(atctb), O(o%) contributions to the quantity Ap in the limit of heavy scalar quarks, i.e. 
we consider the contributions of a Two-Higgs-Doublet model with MSSM restrictions. 
The full analytic result for arbitrary values of the lightest CP-even Higgs boson mass is 
presented. The numerical effects of the leading electroweak MSSM two-loop corrections 
on the precision observables M\y and sin 2 8 e g are analyzed. The electroweak two-loop 
contribution to Myy amounts up to —12 MeV and up to +6 x 10~ 5 for sin 2 # e fj. The 
corrections from the bottom quark loops can become important for large values of 
tan/3. They enter with a different sign than the C(a 2 ) corrections. We furthermore 
investigate the current sensitivity of the electroweak precision observables to the top 
Yukawa coupling in the SM and the MSSM. The prospects for indirectly determining 
this coupling at the next generation of colliders are discussed. 
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1 Introduction 



Theories based on Supersymmetry (SUSY) U are widely considered as the theoretically 
most appealing extension of the Standard Model (SM). They are consistent with the ap- 
proximate unification of the gauge coupling constants at the GUT scale and provide a way 
to cancel the quadratic divergences in the Higgs sector hence stabilizing the huge hierarchy 
between the GUT and the Fermi scales. Furthermore, in SUSY theories the breaking of the 
electroweak symmetry is naturally induced at the Fermi scale, and the lightest supersym- 
metric particle can be neutral, weakly interacting and absolutely stable, providing therefore 
a natural solution for the Dark Matter problem. 

Supersymmetry predicts the existence of scalar partners fi, fn to each SM chiral fermion, 
and spin-1/2 partners to the gauge bosons and to the scalar Higgs bosons. So far, the 
direct search for SUSY particles has not been successful. One can only set lower bounds of 
(9(100) GeV on their masses 0. Furthermore, contrary to the SM two Higgs doublets are 
required resulting in five physical Higgs bosons j3| . The direct search resulted in lower limits 
of about 90 GeV for the neutral Higgs bosons and about 80 GeV for the charged ones ||. 

An alternative way to probe SUSY is via the virtual effects of the additional particles 
to precision observables. This requires a very high precision of the experimental results as 
well as of the theoretical predictions. The most prominent role in this respect plays the 
p-parameter ||. The radiative corrections from vector boson self-energies to the quantity 
Ap constitute the leading, process independent corrections to many electroweak precision 
observables, such as the prediction for Ar, i.e. the M w — M z interdependence, and the 
effective leptonic weak mixing angle, sin 2 6* eff . 

The radiative corrections to the electroweak precision observables within the Minimal 
Supersymmetric Standard Model (MSSM) stemming from scalar fermions, charginos, neu- 
tralinos and Higgs bosons have been discussed at the one-loop level in Refs. ||[7fl, providing 
the full one-loop corrections. More recently also the leading two-loop corrections in 0(aa s ) 
to the quark and scalar quark loops for Ap have been obtained f| as well as the gluonic 
two-loop corrections to the Mw — Mz interdependence 0. Contrary to the SM case, these 
two-loop corrections turned out to increase the one-loop contributions, leading to an en- 
hancement of up to 35% p. 

In this paper we present the leading two-loop corrections to Ap at 0(af ), 0(a t ab), 0(al), 
i.e. the leading two-loop contributions involving the top and bottom Yukawa couplings. 
These contributions are of particular interest, since they involve corrections proportional to 
m\ and bottom loop corrections enhanced by tan f3, the ratio of the two vacuum expectation 
values. For a large SUSY scale, Msusy 3> Mz, the contributions from loops of SUSY 
particles decouple from physical observables fl0|,j8ll. Therefore, focusing on the case of large 
^susy, we derive the leading electroweak two-loop corrections in the limit where besides the 
SM particles only the two Higgs doublets of the MSSM are active. 

As a first step, in Ref. |TTJ we have calculated the 0(a^) corrections in the limit where 



the lightest CP-even Higgs boson mass vanishes, i.e. rrih — > 0. The numerical effect of these 
corrections turned out to be relatively small. However, for the corresponding SM result it 
was found that the Mj| M = limit is only a poor approximation of the result with arbitrary 
M^ M [12]. Since a similar behavior can be expected for the MSSM, we perform the calculation 



of the leading electroweak two-loop corrections, O(a^), 0(a t ab), and O(a^), for arbitrary 
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values of m^. The result obtained in the MSSM is compared with the corresponding SM 
correction of 0(a 2 ) ||12|| . The resulting shift in Myy and sin 2 9 e g is analyzed numerically. 

Since the top Yukawa coupling enters the predictions for the electroweak precision ob- 
servables at lowest order in the perturbative expansion at 0(a 2 ), these contributions allow to 
study the sensitivity of the precision observables on this coupling. Using a simple approach 
in which we treat the top Yukawa coupling in the SM and the MSSM as a free parameter, 
we study the current sensitivities of the electroweak precision observables as well as the 
prospective accuracies at the next generation of colliders. 

The rest of the paper is organized as follows: in Sect. 2 we review the SM and MSSM 
corrections to the quantity Ap and present the details of the calculation of the 0(a 2 ), 
0(a t otb), @( a b) corrections. Explicit formulas for the results of 0(a 2 ), 0(a t ab), and C(a 2 ) 
can be found in Sect. 3 and the appendix. The numerical analysis is performed in Sect. 4. In 
Sect. 5 we analyze the sensitivity of the electroweak precision observables to the top Yukawa 
coupling. We conclude with Sect. 6. 



2 Calculation of the 0(ql 2 ), 0(QL t QL h ), and 0{oL 2 h ) 
corrections 

2.1 One- loop results 

The quantity Ap, 

parameterizes the leading universal corrections to the electroweak precision observables in- 
duced by the mass splitting between fields in an isospin doublet ||. S^ iH /(0) denote the 
transverse parts of the unrenormalized Z and W boson self-energies at zero momentum 
transfer, respectively. Ap gives the dominant contribution to electroweak precision observ- 
ables, such as the W boson mass, Mw, and the effective leptonic mixing angle, sin 2 9 c s. The 
induced shifts are in leading order given by (with 1 — s 2 ^ = c 2 ^ = M^/M§) 

„2 „2 



Mw C W a„ X „;„2 a ^ C W S W 



5M W w 7 2 w 2 Ap, 5 sin 2 8 eS w - 2 w w 2 Ap. (2) 



-w w °w 

In the SM the dominant contribution to Ap at the one-loop level is given by the t/b 
doublet due to their large mass splitting. It reads 

Apr = ^- 2 F (mlml), (3) 



with 



2 x y x 

F (x,y)=x + y log-. (4) 

x-y y 

F Q has the properties F (m 2 ,m 2 ) = F (m 2 ,m 2 ), F (m 2 ,m 2 ) = 0, F (m 2 ,0) = m 2 . Therefore 
for m t ^> rrib eq. (|3|) reduces to the well known quadratic correction in m t , 

s M = _3Gf_ 2 _ 
H0 8V2n 2 * V ' 
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Within the MSSM the dominant correction from SUSY particles at the one-loop level 
arises from the scalar top and bottom contribution to eq. ([]]). For rrif, 7^ it is given by 

Ap " USY = ^Tfe I " sin2 ^ cos2 W K> ™U - sin2 °i cos2 WH' m V 

+ cos 2 6~ t cos 2 e- h F {m 2 k ,m- bi ) + cos 2 6- t sin 2 ^F (m| , m~ 2 ) 
+ sin 2 0- t cos 2 9- b F (m| , mf J + sin 2 £ sin 2 %F (m| , m? 2 ) ] . (6) 

Here mi m~ h .{i = 1,2) denote the stop and sbottom masses, whereas 9f,6i are the mixing 
angles in the stop and in the sbottom sector. 

2.2 Results beyond the one-loop level 

Within the SM the one-loop 0(a) result has been extended in several ways. The dominant 
two- loop corrections at 0(aa s ) are given by |~3 



ApSM^ = -Ap^- (1 + ^/3). (7) 

6 7T 

These corrections screen the one-loop result by approximately 10%. Also the three-loop 
result at 0(aa 2 ) is known. Numerically it reads |L4| 



Apf = -^^m 2 f^) 2 • 14.594... . (8) 

Furthermore the leading electroweak two-loop top quark contributions of 0(a 2 ) have been 
calculated. They enter the electroweak precision observables together with the one-loop 
contribution according to 

P= I^V A P = A Po + A Pi- (9) 



First the result for Api in the limit Mj| = had been evaluated [15 



6!\£f H =o = 19-27T 2 . (10) 



Later the full 0(a 2 ) result without restrictions in the Higgs boson mass became available fi~2|l . 
where Sf^ H=0 extends to 

S^ h¥0 = 19 - 2tt 2 + fct(m t , M H ) . (11) 

Here ict(mt, Mh) contains the extra terms arising from a non-vanishing Higgs boson mass. 
Recently also first electroweak three- loop results in the limit of Mh = became available ||16|| . 
Numerically they read 

SM.af _ ( Gp ^ 3 



AP 2 ,ZU = \ • 2J,I7J - (12) 



A «S - (j^f) 2-9394. (!3) 

In the MSSM up to now the two-loop calculations have been restricted to the leading 
0(aa s ) corrections to the scalar quark loops ||. They consist of the rather lengthy result 
for gluino exchange, which decouples for trig — > oo, and of the compact correction for the 
gluon exchange contribution [RJ: 



A^SUSY _ G F 



' gluon 4V27T 2 7r 



- sin 2 % cos 2 9- t Fi(m\ , m| ) - sin 2 6 b cos 2 0~ b Fi (m^ , m~ 2 ) 



+ cos 6* t -cos 9- b Fi{m^ m^ ) + cos ^sin ^Fi(m^,m 



+ sin 2 t - cos 2 0gi5\ (m| , ) + sin 2 t - sin 2 (m| , m| ) , (14) 



with 



F^x.y) = x + y-2 log - 



2 + - In - 

2/ 2/ 



(x + y)x 2 , nX „. , T . / x\ 

+ 7 lo S - - 2(x - 2/)Li 2 1 - - , (15) 

[x-yY y V Vj 

where F\ has the properties Fi(m 2 ,m 2 ) = Fi(m 2 ,m 2 ), Fi(m 2 ,m 2 ) = 0, i<i(m 2 ,0) = m 2 (l + 
vr 2 /3). 

Contrary to the SM case where the strong two- loop corrections screen the one- loop result, 
the 0(aa s ) corrections in the MSSM increase the one-loop contributions by up to 35%, thus 
enhancing the sensitivity to scalar quark effects. Another difference between the SM and 
the MSSM is the m t dependence of the leading contribution to Ap within the SM. They are 
~ m 2 for the one-loop and ~ m\ for the two-loop correction leading to sizable shifts to the 
precision observables. Concerning the corrections from loops of SUSY particles, on the other 
hand, no large prefactor ~ m 2 at the one- loop level is present. This behavior changes with 
the leading electroweak two-loop SUSY corrections, which are ~ m|, i.e. of 0(a 2 ). Therefore 
we concentrate on these and the corresponding 0(a t ab), 0{a b ) corrections in this paper. 
Since the SUSY loop contributions in the MSSM decouple if the general soft SUSY-breaking 



scale goes to infinity, Msusy — > oo [10|,|8|, the leading contributions for large Msusy arise 



from a Two-Higgs-Doublet model with MSSM restrictions. 
2.3 The Higgs sector of the MSSM 

Contrary to the SM, in the MSSM two Higgs doublets are required [|J. At the tree-level, 
the Higgs sector can be described with the help of two independent parameters (besides g 
and g'): the ratio of the two vacuum expectation values, tan/? = v 2 /vi, and Ma, the mass 
of the CP-odd A boson. The diagonalization of the bilinear part of the Higgs potential, i.e. 
the Higgs mass matrices, is performed via orthogonal transformations with the angle a for 
the CP-even part and with the angle /3 for the CP-odd and the charged part. The mixing 
angle a is determined through 

, M'\ • M i 7T 

tan2a = tan2/? M ^_ M f ; -- < a < . (16) 
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One gets the following Higgs spectrum: 



2 neutral bosons, CV = +1 
1 neutral boson, CV = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



h°,H° 
A 

H + ,H- 

s~iO r~i— 



(17) 



At the tree level, the Higgs boson masses expressed through Mz, Mw and Ma are given by 

(18) 



m h 


1 

2 


Mi + Mf- V 


l{M\ + M\Y 


- AM\M 2 Z cos 2 2(3 




1 

2 


M\ + M 2 Z + v 


/(Mi + M 2 ) 2 


- 4M|M| cos 2 2/3 


m 2 H± 


= M 


\ + M 2 w 






m G 


= Ml 






m 2 G± 


= M*,, 







(19) 

(20) 
(21) 
(22) 

where the last two relations, which assign mass parameters to the unphysical scalars G° and 
are to be understood in the Feynman gauge. 



2.4 Evaluation of the (D(cq), 0(a t ab), 0(al) contributions 

In order to calculate the 0(a 2 ), 0(actatb), 0( a 1) corrections to Ap, see eq. (|T|), the Feynman 
diagrams generically depicted in Fig. have to be evaluated for the Z boson (V = Z) and the 
W boson (V = W) self-energy. We have taken into account all possible diagrams involving 
the t/b doublet and the full Higgs sector of the MSSM, see Sect. |2~3| . 

The two-loop diagrams shown in Fig. [TJ have to be supplemented with the corresponding 
one-loop diagrams with subloop renormalization, depicted generically in Fig. [| The coun- 
terterms that enter the calculation are the top mass counter term, 5mt, the Higgs boson 
mass counter term, <5M|, and the tadpole counter terms, ST^ and 5Th- The renormalization 
constants have been derived in the on-shell scheme as outlined in Ref. [0. The wave func- 
tion renormalization constants, entering via the diagrams in Fig. 0, drop out as required. 
The Feynman diagrams for the insertions of the fermion and Higgs mass counter terms are 
shown in Fig. [| 

The amplitudes of all Feynman diagrams, shown in Figs. [l|-|3], have been created with the 
program FeynArtsS [|18|| , making use of the MSSM model file |T9[] (where only the non-SM like 



counter terms had to be added). The reduction to scalar integrals has been performed with 
the program TwoCalc, based on the reduction method of Ref. |2(|. As a result we obtained 



the analytical expression for Ap depending on the one-loop functions Aq and Bq |2T| and on 
the two-loop function T 134 |20|,^2j. For the further evaluation the analytical expressions for 
A , B and T 134 have been inserted. 

In order to derive the leading contributions of 0(a 2 ), 0(a t aib) an d 0(a 2 ) we extracted 
the contributions proportional to y 2 , y t yb and y 2 , where 

V2m t V2m b 

Vt = — '■ 7T> Vb = w ■ ( 23 ) 

v sin (3 v cos [3 
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q q 




Figure 1: Generic Feynman diagrams for the vector boson self-energies, 
{V = {Z, W}, q = {t, b}, 0, X = {h, H, A, H±, G, G±}). 




q 



Figure 2: Generic Feynman diagrams for the vector boson self-energies with counter term 
insertion, (V = {Z, W}, q = {t, b}, <j>, X = {h, H, A, H ± , G, G ± }). 



q t) 




q q 



Figure 3: Generic Feynman diagrams for the counter term insertions, 
(q = {t,b}, ( j ) ={h,H,A,H ± ,G,G ± }). 
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The coefficients of these terms could then be evaluated in the gauge-less limit, i.e. for 
M w , M z -> (keeping c w = M w /M z fixed). 

For the Higgs masses appearing in the two-loop diagrams we use the following relations, 
arising from the gauge-less limit 

m 2 H± = M\ , 
m 2 G = , 

m 2 G± = . (24) 

Applying the corresponding limit also in the neutral CP-even Higgs sector would yield for 
the lightest CP-even Higgs-boson mass ml = (and furthermore m 2 H = M\, sin a = 
— cos (3 , cos a = sin/5 ). Since within the SM the limit Mf^ 1 —>■ turned out to be only 
a poor approximation of the result for arbitrary M^ M , we keep in our calculation a non- 
zero m\ (which formally is a higher-order effect). Keeping m^ as a free parameter is also 
relevant in view of the fact that the lightest MSSM Higgs boson receives large higher order 
corrections f2"5fl , which shift its upper bound up to 135 GeV (for M SU sy < 1 TeV and 
m t = 175 GeV) []17|,|24|. These corrections can easily be taken into account in this way 



(in the Higgs contributions at one-loop order, however, the tree-level value of mh should be 
used). Keeping a arbitrary is necessary in order to incorporate non SM-like couplings of the 
lightest CP-even Higgs boson to fermions and gauge bosons. 

On the other hand, keeping all Higgs-sector parameters completely arbitrary is not pos- 
sible, as the underlying symmetry of the MSSM Lagrangian has to be exploited in order 
to ensure the UV-finiteness of the two-loop corrections to Ap. We thus have enforced only 
those symmetry relations in the neutral CP-even Higgs sector which are explicitly needed in 
order to obtain a complete cancellation of the UV-divergences. 

In the following, we separately consider the 0(a 2 ) corrections, corresponding to the limit 
where yb = 0, and the full 0(a 2 ), 0(a t ab), 0(a 2 ) contributions. The 0(a 2 ) corrections are 
by far the dominant subset within the SM, i.e. the 0(a t ab) and 0(a 2 ) corrections can safely 
be neglected within the SM. The same is true within the MSSM for not too large values 
of tan/3. Thus, we first consider the scenario where only the 0(a 2 ) corrections need to be 
taken into account and then discuss the result in the case where the 0{a.tOib) and 0(a 2 ) 
corrections are non- negligible. 

In the case of the 0(a 2 ) corrections, no further relations in the neutral CP-even Higgs 
sector are necessary, i.e. we keep the parameters mh, run and a arbitrary in the evaluation 
of the 0(a 2 ) corrections. For these contributions also the top Yukawa coupling y t can be 
treated as a free parameter, i.e. it is not necessary to use eq. (^). As a consistency check of 
our method we recalculated the corresponding SM result jTT2| and found perfect agreement. 

Concerning the corrections to Ap with yb ^ 0, the SU(2) structure of the fermion doublet 
requires further symmetry relations. Within the Higgs boson sector it is necessary, besides 
using eq. (|24|), also to use the relations for the heavy CP-even Higgs boson mass and the 
Higgs mixing angle, 

m 2 H = M 2 A , 
sin a = — cos (3 , 

cos a = sin/3 . (25) 
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On the other hand, can be kept as a free parameter. The couplings of the lightest 
CP-even Higgs boson to gauge bosons and SM fermions, however, become SM-like, once 
the mixing angle relations, eq. (|25|), are used. Furthermore, the Yukawa couplings can no 
longer be treated as free parameters, i.e. eq. fl23|) has to be employed, which ensures that 
the Higgs mechanism governs the Yukawa couplings. Also in this case we have evaluated 
the corresponding SM corrections. As expected, the corrections arising from y^ ^ are 
numerically insignificant within the SM. 



The analytical result for the 0(oq) contributions for 
the special case = 



For illustration, we discuss in this section the result for the 0(af) contributions for the 
special case where = 0. The result for the 0(af) corrections for arbitrary parameters in 
the CP-even Higgs sector can be found in Appendix A.l, while the result for the full 0(af), 
0(a t otb) and 0(af) corrections can be found in Appendix A. 2. The full results have been 
included into the code FeynHiggs |fj5 . 



In order to simplify the expression for the 0(af) contributions as far as possible, we use 
in this example the relations eqs. (^4j) and (^) as well as rrih = 0. The only remaining 
parameters in this case are the top quark mass, m t , the CP-odd Higgs boson mass, M4, and 
tan (3 (or sp = tan 



The analytical result obtained as described in Sect. can conveniently be expressed in 
terms of 



A = 



(26) 



The two-loop contribution to the p -parameter then reads: 



Ap* 



SUSY 
Higgs,m h =0 



ri2 

4 

— 7 "V 

A 1 



128 n 4 



A 2 



x 



|li 2 ((l — VI -4 A) /2 



-2Li 2 U-t) [5-14A + 6A 2 ] 



=A 



+ \og 2 (A) 



1 + 



y/l-AA 



-.A 



log 2 



1 - y/l-AA 



VI - 4A 



log(A) 



=A 



2 - 20 A 



io ; , ( L-^O?) vr^4A(i -2 a) 



.1 + VI 
-log(|l/A-l|) {A -If 



2 ^™A + l --2A^ 
-3 + 12 A 3 1-s 2 



- 11 A + 19 : 



A 2 



1 - s 



2 ( ' 



(27) 
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with 

A = 3 - 13A+ 11 A 2 . 
In the limit of large tan/3 (i.e. (1 — s|) <C 1) one obtains 



:28i 



a SUSY 



m h =0 



r 2 

3 — — — r mi 



128 7T 



4 * 



19 



27T 2 + 0(l-4) 



/3 



(29) 



Thus for large tan /5 the SM limit with Mjj — > is reached. 

In order to investigate the decoupling behavior of Ap^^ s mh=0 , the result for Apf^ gsmh=0 
in eq. ( ^7f) can be expanded for small values of A, i.e. for large values of Ma- 



a sus 

^ri,Hi 



SUSY 

ggs,m h =0 



4 

r 4 * 



128 7T 4 



19 - 2n z 



1-4 



log 2 A + y ) (&4 + 32A 2 + 132A 3 + 532A 4 ) 



+ log(A)— (560A + 2825 A 2 + 11394A 3 + 45072A 4 ) 

oU 

— (2800A + 66025A 2 + 300438A 3 + 1265984A 4 ) + O (A 5 ) 

1800 V ) \ ) 



In the limit A — » one obtains 



SUSY 
,Higgs,m^=0 



128 7T 4 



[19-2 7T 2 ] 



(30) 



(31) 



i.e. exactly the SM limit for M^ M — > is reached. This constitutes a consistency check, since 
in the limit A — > the heavy Higgs bosons are decoupled from the theory. Thus only the 
lightest CP-even Higgs boson should remain, which in the 0(a 2 ) approximation (neglecting 
higher-order corrections) has the mass rrih = 0. 

An expansion for small values of Ma as well as an analysis of the quality of these expan- 



sions can be found in Ref. ITT 



The more general expressions, i.e. with rrih 7^ 0, at 0(af), 0(a t ab), and C(a 2 ) can be 
found in the appendix. 



4 Numerical analysis 

4.1 The 0{af) contributions 
4.1.1 Comparison for Ap 

In Fig. [| the size of the leading 0(a 2 ) MSSM corrections, eq. fl3"7|), is compared with the 
leading 0(af) contribution in the SM (jTJ|, with the leading MSSM corrections arising from 
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the i/b sector at 0(a) |7j, and with the corresponding gluon-exchange contributions of 
0(aa s ) || (the 0(aa s ) gluino-exchange contributions |§, which go to zero for large rrig, have 
been omitted here). The numerical effects of the different contributions to Ap are shown as 
a function of a common SUSY mass scale, Msusy (which enters the diagonal entries in the 
t mass matrix). For the MSSM parameters we have chosen the values as specified in the 
m max bendrmark scenario [26], i.e. X t = 2 Msusy ; where m t X t is the off-diagonal entry in the 
t mass matrix. For our conventions in the t sector, see Ref. [17]. The other parameters are 



p = 200 GeV, A b = A t . p is the Higgs mixing parameter and A t)b are the trilinear Higgs-i, b 
couplings, respectively. While Ma has been set to Ma = 300 GeV, for tan (3 we have chosen 
two typical values, tan/3 = 3 as a low and tan/3 = 40 as a high value. (Smaller tan/3 values 
within the m™ ax scenario, where m t is fixed to m t = 174.3 GeV and Msusy < 1000 GeV, 
are disfavored by the LEP Higgs boson searches P),|2"7fl.) From these parameters the values 
for rrih, triH and a have been obtained. For the numerical evaluation of the CP-even Higgs 



boson sector, we have used the results from the t/t sector as presented in Refs. fT7|,[2"8, 
The SM O(o.t) corrections, eq. (0), have been evaluated using the result of rrih as the SM 
Higgs boson mass. Msusy enters the MSSM 0(af) corrections (where as described above, 
we have neglected the SUSY loop contributions) only indirectly via its effect on m h . 

Fig. |] shows the decoupling of the effects of scalar quark loops with increasing Msusy- 
The 0(aa s ) SUSY corrections are always about an order of magnitude smaller than the 0(a) 
squark loop contributions to Ap. The decoupling with Msusy indicates that for large values 
of Msusy the contributions from quarks and scalar quarks within the MSSM essentially 
reduce to the quark loop corrections. This motivates to approximate the 0(af) corrections 
in the full MSSM by the Two-Higgs-Doublet model part. 

The 0{a\ ) corrections involving quarks and the Higgs sector of the MSSM turn out to be 
larger than the 0(aa s ) SUSY corrections for all values of Msusy ~ 200 GeV. This is related 
to the enhancement by the prefactor mf/Myy- The 0(af) corrections even exceed the 0(a) 
squark loop corrections for M S usy ~ 600 GeV, i.e. these contributions can compensate each 
other as they enter with different sign. This applies also for the no-mixing scenario (X t = 0, 
Msusy = 2000 GeV) |2§, which is not shown here. 

In Fig. ^ we analyze the dependence of the 0(af) contributions to Ap on the lightest 
CP-even Higgs boson mass, rrih. For the MSSM parameters we have again chosen values as 
specified for the m™ ax and the no-mixing scenario. While tan/3 has been fixed to tan/3 = 
3,40, the CP-odd Higgs boson mass has been varied from 50 GeV to 1000 GeV. 

As can be seen in Fig. ||, the 0(af ) MSSM contribution is of O(10~ 4 ). It is always larger 
than the corresponding SM result. In the limit of large Ma, i.e. at the endpoint of the rrih 
spectrum, the difference of the SM and the MSSM result are numerically very small. This 
is in accordance with the decoupling behavior that we have discussed analytically for the 
special case with rrih = 0, see Sect. || 

In Fig. |6| the decoupling behavior for large Ma of Ap^^ is analyzed numerically. Since 
Ma is the only non-SM scale that directly enters Ap^^ s , the result should become SM- 
like once Ma is taken to very large values. This is shown in the m™ ax and the no-mixing 
scenario for tan/3 fixed to tan/3 = 3(40) in the left (right) plot of Fig. [| The SM value 
of Ap™ iggs is obtained for M^ M = m h . While for the small tan/? value the decoupling is 
very slow and is reached only for Ma ^ 3 TeV, for the large tan (3 value Ap^^ decouples 
already for very small Ma values. This numerical behavior is analogous to the analytical 
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Figure 4: The contribution of the leading 0(a%) MSSM corrections, Apf^ gs , eq. fl37D , 
is shown as a function of Msusy for Ma = 300 GeV and tan/5 = 3 (left plot) or tan/3 = 
40 (right plot) in the m™ ax scenario. Apf £ is compared with the leading 0(ot£) SM 
contribution and with the leading MSSM corrections originating from the t/b sector of 0(a) 
and 0(aa s ). Both O{o%) contributions are negative and are for comparison shown with 
reversed sign. In the right plot the O(o£) corrections differ by about 1.5 x 10~ 7 , which is 
not visible in the plot. 





110 120 

m h [GeV] 




110 1 20 

m h [GeV] 



Figure 5: The contribution of the leading 0(a^) MSSM corrections, Apf^ gs , is shown as 
a function of rrih for tan (3 = 3 (left plot) and for tan (3 = 40 (right plot) in the m™ ax and the 
no-mixing scenario. Apf is compared with the leading 0(af) SM contribution. 
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Figure 6: The 0(af) MSSM contribution to Ap in the m™ ax and the no-mixing scenario 
is compared with the corresponding SM result with M^ M = m^. is obtained in the left 
(right) plot from varying Ma from 100 GeV to 3000(300) GeV, while keeping tan/3 fixed at 
tan/5 = 3(40). 

result described in Sect. |3| for the — > limit. However, it should also be noted that for 
small values of Ma the behavior of the SM and the MSSM contributions is very different. 
While the SM contribution depends sensitively on M| M , in the MSSM for M A > 100 GeV 
(corresponding to nih > 90 GeV, see Fig. H) the dependence on the Higgs boson masses is 
much less pronounced, see in particular the right plot of Fig. ||. 

4.1.2 Effects on precision observables 

In this section the numerical effect of the 0(af) corrections on the electroweak precision 
observables, My/ and sin 2 9 e s is analyzed. In addition to the MSSM 0(af ) correction to 
5My/ and 5 sin 2 # e fj, we also present the effective change from the SM result (where the SM 
Higgs boson mass has been set to m^) to our new MSSM result. 

In Fig. [7] the absolute contribution and the effective change for the W boson mass is 
presented. For the numerical evaluation we have chosen the m™ ax benchmark scenario (where 
A^susy = 1000 GeV). In the left plot tan/5 is fixed to tan/3 = 3,40, while Ma is varied from 
50 GeV to 1000 GeV, resulting in the Higgs boson mass m h . The effect of the 0(a 2 t ) MSSM 
contributions on 8My/ amounts up to —12 MeV. For large Ma and/or large tan f3 it saturates 
at about —10 MeV. The effective change in My/ is significantly smaller. It amounts up to 
—3 MeV and goes to zero for large Ma as expected from the decoupling behavior. In the right 
plot of Fig. [7| 5My/ is shown as a function of tan f3. Ma is kept fixed to Ma = 100, 300 GeV. 
The effect of Apf^^ saturates for large tan (3. For a small CP-odd Higgs boson mass, 
Ma = 100 GeV, a shift of —2 MeV in My/ remains also in the limit of large tan (3, since the 
two Higgs doublet sector does not decouple from the MSSM. For large Ma, Ma = 300 GeV, 
for nearly all tan (3 values the effective change in My/ is small. 



12 



0.002 



0.000 



-0.002 



> -0.004 
O 



-0.006 



-0.008 



-0.010 



-0.012, 



ii.i! 

- 


i i i i I i 








m h™ X 




tanp = 3 _ 




tanp = 40 




tanp = 3 (MSSM - SM) _ 




tanp = 40 (MSSM - SM) - 








' " 



90 



100 



110 120 130 140 

m, [GeV] 



0.002 



0.000 



-0.002 



> -0.004 



-0.006 



-0.008 



-0.010 - 



-0.012 



_' ' ' ' I 




- 
- 






'! 




- 


1 




_ 


- 


m h max 


- 




M fl = 1 00 GeV 






M A - 300 GeV 






M fl -100 GeV (MSSM 


SM) - 




M A = 300 GeV (MSSM 


SM) ! 










i i I i i i i I i i i i I i 





10 



20 30 

tanp 



40 



50 



Figure 7: The absolute O(af) MSSM contribution and the effective change in SM W is shown 
for Msusy — 1000 GeV in the m™ ax scenario. The other parameters are fx = 200 GeV, = 
A t . ffift is obtained in the left (right) plot from varying M\ from 50 GeV to 1000 GeV, while 
keeping tan f3 fixed at tan /3 = 3, 40 (from varying tan /3 from 2 to 40, while keeping Ma fixed 
at M A = 100, 300 GeV.) 
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Figure 8: The absolute 0{af) MSSM contribution and the effective change in <5sin 2 ^ e ff 
is shown for M S usy = 1000 GeV in the m™ 3 ^ scenario. The other parameters are fx = 
200 GeV, Aj, = A t . rrih is obtained in the left (right) plot from varying Ma from 50 GeV to 
1000 GeV, while keeping tan f3 fixed at tan f3 = 3, 40 (from varying tan /3 from 2 to 40, while 
keeping M A fixed at M A = 100, 300 GeV.) 
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The absolute contribution and the effective change for 8 sin 2 9 e ^ is shown in Fig. || for 
the same parameters as in Fig. [7]. The absolute effect is around +6 x 10~ 5 . The effective 
change ranges between +3 x 10~ 5 for small tan/? and small Ma and approximately zero for 
large tan (3 and large Ma- 

The effects of the 0(af) MSSM corrections in My/ and sin 2 9 e s discussed above are smaller 
than the current experimental errors, SM^ P = 34 MeV and Ssxn 9^ = 17 x 10~ 5 ||30|| . 
However, their inclusion is crucial in order to reduce the theoretical uncertainties from 
unknown higher order corrections within the MSSM to a similar level as in the SM of 

±7 x 10 -5 pi],p2J. Achieving this level of theoretical 



5M$ eo Rj ±5 MeV and 5 sin 2 9^ co 



accuracy will be mandatory in particular in view of the prospective accuracies at a future 



linear collider running on the Z peak and the WW threshold (GigaZ), 5M' W 



exp,fut 



7 MeV 



and 5 sin 2 fl e e * p ' fut 



1.3 x 10" 



4.2 The 0(o$), 0(a t at,), and 0(af) contributions 

In this section the numerical effect of the 0(af), 0(a t ab), and 0{al) corrections on Ap 



is analyzed. As discussed in Sect. [2~4| , for these corrections is was necessary to employ the 
Higgs sector restrictions as given in eqs. (|24l) and fl25|). This implies that the couplings of 
the lightest CP-even Higgs boson to gauge bosons and SM fermions are SM-like. Corrections 
enhanced by tan (3 thus arise only from the heavy Higgs bosons, while the contribution from 
the lightest CP-even Higgs boson resembles the SM one. 



0.0000 



-0.0001 



-0.0002 



-0.0003 



90 



tanp = 40 
a (SM) 

a (MSSM, no mix.)) 
a 2 (MSSM, m h max ) 



100 



110 120 

m h [GeV] 



130 



140 



-0.0001 Or 



-0.00012 - 



-0.00014 



M A = 300 GeV 

a (SM, M H = m™ m ") 

a 2 (MSSM, no mix.)) 

a (SM, M H = m h max ) 

a (MSSM, m h max ) 




-0.00018 



-0.00020 I 



10 



20 30 

tanp 



40 



50 



Figure 9: The 0(a 2 ), 0(a t a b ), and £>(a 2 ) MSSM contribution to Ap in the m^ ax and the 
no-mixing scenario is compared with the corresponding SM result with Mj| M = m^. In the 
left plot tan/3 is fixed to tan/? = 40, while Ma is varied from 50 GeV to 1000 GeV. In the 
right plot Ma is set to 300 GeV, while tan (3 is varied. The bottom quark mass is set to 
m b = 4.25 GeV. 

In Fig. ^| we show the result for the C(a 2 ), 0(a t ab), and C(a 2 ) MSSM contributions to 
Ap in the m™ ax and the no-mixing scenario, compared with the corresponding SM result 
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with M^ M = m/j. In the left plot tan/? is fixed to tan/3 = 40 and M4 is varied from 50 GeV 
to 1000 GeV. In the right plot M A is fixed to M A = 300 GeV, while tan (3 is varied. 

For large tan/3 the 0(a t &b) an d 0(al) contributions yield a significant effect from the 
heavy Higgs bosons in the loops, entering with the other sign than the 0{oc%) corrections, 
while the contribution of the lightest Higgs boson is SM-like. As one can see in Fig. |9], for 
large tan/3 the MSSM contribution to Ap is smaller than the SM value. For large values 
of Ma, the SM result is recovered. The effective change in the predictions for the precision 
observables from the 0(a t Qb) an d 0(af) corrections can exceed the one from the 0(a^) 
corrections. It can amount up to 5M W rs +5 MeV and <5sin 2 # eff ~ — 3 x 10 -5 for tan (3 = 40. 

5 Constraints on the top Yukawa coupling in the SM 
and the MSSM 

a 2 

The Ap x jj. corrections in the SM and the MSSM are of particular interest, since these are 
the leading corrections in which the top and bottom Yukawa couplings, i.e. the coupling of 
Higgs bosons to top and bottom quarks, enter the predictions for the electroweak precision 
observables. Thus, the electroweak precision tests of the SM and the MSSM provide some 
sensitivity to the Yukawa couplings in these models. 

In order to exemplify this sensitivity, we use a simple approach in which we treat the top 
Yukawa coupling in the SM and the MSSM as a free parameter. While a complete calculation 
of top and bottom contributions, as discussed in the previous sections, requires the relation 
between the Yukawa coupling and fermion mass within the SM and the MSSM, this relation 
is not formally needed if one restricts to the top contributions only. Numerically, this is an 
excellent approximation within the SM and also in the MSSM for not too large tan j3. 

In the following we analyze the sensitivity to the top Yukawa coupling in the SM and 

2 

the MSSM. Since in the MSSM contributions beyond the Ap"' Hi corrections are not yet 
known, for this comparison we restrict the SM contributions also to the leading electroweak 

2 

^Pi Higgs term [JT2] , neglecting the formally subleading electroweak two-loop corrections to 



the precision observables which can, however, be of similar size. 

Fig. [H] shows the effect of varying the top Yukawa coupling in the SM and the MSSM 
for the precision observables Mw and sin 2 9 e g in comparison with the current experimental 
precision. The allowed 68% and 95% C.L. contours are indicated in the figure. The Yukawa 
coupling is scaled in the following way, 

yt = xy? M J 0<x<3, (32) 

and analogously in the MSSM. A shift of this kind in the relation between a fermion mass and 
the corresponding Yukawa coupling can occur for instance in the MSSM (see e.g. Ref. ||35||), 

Vim 1 /qq\ 

Vt = — '■ 77" : , a > 33 
v sin /3 1 + A t 

where A t is induced by SUSY loop corrections. Here we do not assume any particular scenario 
but use the variation of the top Yukawa coupling only for demonstrating the sensitivity to 
this parameter. 
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Figure 10: The effect of scaling the top Yukawa coupling in the SM (upper plot) and the 
MSSM (lower plot) for the precision observables M w and sin 2 # eff is shown in comparison with 
the current experimental precision. The variation with m t and Aa^ad is shown within their 
current experimental errors. For the SM evaluation, M| M has been set to the conservative 
value of M^ M = 114 GeV (see text). For the MSSM evaluation the parameters are Msusy = 
1000 GeV, X t = 2000 GeV, M A = 175 GeV, tan (3 = 3 and fi = 200 GeV, resulting in 
m h w 114 GeV. The obtained upper bounds are y t < 1.3 yf M and y t < 1.7^ MSSM for 
m t = 174.3 GeV as well as y t < 2.2 yf M and y t < 2.5^ MSSM for m t = 179 A GeV, all at the 
95% CL. 
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For the evaluation of M^y and sin 2 9 e s in the SM and the MSSM, we take into account 
the complete one-loop results as well as the leading two-loop 0(aa s ) and 0(af) corrections 

oil 



fas discussed in Refs. 



Since the SM prediction deviates more from the experimental 
central value for increasing values of Mjj , we have chosen in the figure Mjy =114 GeV |[7|] 
as a conservative value. The current la uncertainties in mt and Acthad are also taken into 
account, as indicated in the figures. Varying the SM top Yukawa coupling (upper plot) yields 
an upper bound of y t < 1.3 yf M for m t = 174.3 GeV and of y t < 2.2 yf M for m t = 179.4 GeV, 
both at the 95% C.L. These relatively strong bounds are of course related to the fact that 
the theory prediction in the SM shows some deviation from the current experimental central 
value. 



The lower plot of Fig. [10] shows the analogous analysis in the MSSM for one particular 
example of SUSY parameters. We have chosen a large value of Msusy, ^susy = 1000 GeV, 
in order to justify the approximation of neglecting the 0(af) contributions from SUSY loops. 
The other parameters are X t = 2000 GeV, M A = 175 GeV, tan (3 = 3 and /i = 200 GeV, 
resulting in rrih ~ 114 GeV (for comparison with the SM case). The SUSY contributions to 
Mw and sin 2 9 e s lead to a somewhat better agreement between the theory prediction and 
experiment and consequently to somewhat weaker bounds on y t . In this example we find 
y t < 1.7?/ t MSSM for m t = 174.3 GeV and y t < 2.5% MSSM for m t = 179.4 GeV, both at the 
95% C.L. 

In order to demonstrate the sensitivity of future colliders for the determination of the 
top Yukawa coupling from electroweak precision observables, we list in Tab. [1| the bounds 
on y t obtainable at the LHC and a future LC with GigaZ option |34| . Here we assume that 
the future experimental central values of Mw and sin 2 9 e g agree with the theory predictions 
for yt/yf M = 1 and yt/yf 18 ^ = 1, respectively. An accuracy in the indirect determination 
of yt of about 40% can be achieved with the GigaZ precision at the 95% C.L. This is similar 
to the accuracy achievable from the ti threshold measurements, see Ref. [p9[| . The results 
in Tab. |T| are the same for the SM and our SUSY example, since the only difference (after 
assuming that the future experimental central values of My/ and sin 2 9 e g agree perfectly with 
the SM or MSSM predictions) are the relatively small deviations at 0(a^) between the SM 
and the MSSM shown in Figs. % |. 
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Table 1: 

5m t = 2 GeV or 1 GeV, 5M W 

we use 5m t = 0.1 GeV, 5M W = 7 MeV and 5 sin 2 8 eS = 1.3 x 10 

assume a future uncertainty of tfActhad = 5 x 10~ 5 . The bounds on yt/yf M and yt/yt 1 ^ 1 ^ are 
given at the 95% C.L., assuming that the theory predictions agree with the experimental 
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^2fl. For A«had we 



central values for M w and sin 9 e s for y t /y- 



SM,SUSY 
t 



1. The SUSY parameters are chosen 



as for Fig. 10 
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6 Conclusions 



We have calculated the leading 0(a 2 ), 0(a t at,), and 0(af) corrections to Ap in the MSSM 
in the limit of heavy squarks. The analytical results for arbitrary values of the lightest CP- 
even Higgs boson mass have been presented. While for the O(a^) result all parameters in the 
CP-even Higgs sector could be kept arbitrary, for the full result further tree-level relations 
had to be employed, which lead to SM-like couplings of the lightest CP-even Higgs boson to 
gauge bosons and SM fermions. 

Numerically we compared the effect of the new MSSM contribution with the leading 
0{o%) SM contribution. For small tan/5, it is sufficient to restrict to the 0(af) corrections. 
Their numerical effect is always larger than the 0{af) SM contribution. The corrections 
to the precision observables Mw and sin 2 # e g- amount up to —12 MeV for Mw and about 
+6 x 10~ 5 for sin 2 # c flf. The effective change from the SM 0(ptf) result with M^ M = is 
smaller. It amounts up to —3 MeV for M w and +2 x 1CT 5 for sin 2 9 cS . This effective change 
goes to zero for large M A , i.e. the non-SM contribution decouples. 

The 0(atab) and 0(al) contributions become important for large tan/3. They enter 
with a different sign than the 0{ot%) corrections and can overcompensate the latter. For 
large tan f3 the effective change in the predictions for the precision observables from the 
whole O(af), 0(a t ab), and O(al) corrections can amount up to 5M W rs +5 MeV and 
5 sin 2 6 eS ~ — 3 x 1CT 5 . For large Ma also in this case the SM result is recovered. 

The MSSM corrections to the electroweak precision observables discussed here are impor- 
tant in order to reduce the theoretical uncertainties from unknown higher order corrections 
within the MSSM to a similar level as currently reached for the SM. Achieving this level of 
theoretical accuracy will be mandatory in particular in view of the prospective accuracies at 
a future linear collider running on the Z peak and the WW threshold. 

We have furthermore discussed the sensitivity of the electroweak precision observables to 
the top Yukawa coupling, which enters at the two-loop level. Varying the SM top Yukawa 
coupling and requiring consistency with the present experimental values of Mw and sin 2 9 e g 
at the 95% C.L. yields an upper bound of y t < 1.3yf M for m t = 174.3 GeV. This bound can 
be relaxed within the MSSM, where additional contributions from SUSY loops to the elec- 
troweak precision observables can lead to a better agreement with the experimental data. 
We have also analyzed the sensitivity of future colliders for the determination of the top 
Yukawa coupling from electroweak precision observables, assuming that the future experi- 
mental central values of Mw and sin 2 9 e g agree with the theory predictions for unmodified 
Yukawa couplings. An accuracy in the indirect determination of y t of about 40% can be 
achieved with GigaZ precision at the 95% C.L., which is similar to the accuracy achievable 
from tt threshold measurements. 
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Appendix 



A Analytical results or arbitrary values of rrih 

A.l The result for the O(a^) contributions for arbitrary parame- 
ters in the CP-even Higgs sector 

We give here the analytical result for the 0(a 2 ) contributions where irih, m# and a are 
kept as arbitrary parameters, see Sect. [2.4| . The full result for the 0(a 2 ) contributions can 
conveniently be expressed in terms of 



A = 



m t 
Ml' 



H = 



m 



2 ' 



mi 



(34) 



"H h 

We furthermore use the abbreviations s x = sin x, c x = cos x and 

3A 2 (1 - 6x + 10x 2 ) - AAx(l - hx + 7x 2 ) + x 2 (l - Ax + 6x 2 ) 
A 2 (-l + 12x + Ylx 2 ) - 2Ax{-\ + lOx + Ax 2 ) + x 2 (-l + 8x) 
2x(-l + Ax) [A(l-Ax) -x] 

2A 2 (-1 + AA)(h - H) \2A 3 + A 2 (-l + 2h + 2H) - 6AhH + hH 
2A(-1 + AA)(h - H) [20A 4 - 2A 3 (5 + 6h + QH) 
+A 2 (1 + 8H + 8h + AhH) — A(h + H + 6hH) + hH 
-2A(-1 + AA)(h - H) [l8A 4 - A 3 (9 + Uh + 1AH) 

+A 2 (1 + 8H + 8h + WhH) — A(h + H + 7hH) + hH 

AA 5 + A 4 (5A + 1360 - 2A 3 (17 + 108/t + 5Ah 2 ) 
+A 2 (5 + 9Ah + U6h 2 ) - 2Ah(6 + 29h) + 7h 2 . 

The two-loop contribution to the p -parameter then reads: 



A, 

a; 

Y x 
Pi 



P 



P, = 



(35) 



a sus 

^Pl,Hi 



SUSY 
ggs 



G 2 F 4 1 
128 7T 4 1 s 2 p 



(36) 



Li 2 ((1-Vl- 4i7)/2 



s 2 a A H + S/3S a (sps a + cpc a )Y H 



Vl-AH(A-H) 2 H 2 

[spCf3(spCfS + S a C a ) + Sf3S a (spS a + CpC a ) + 2s 3 3 S a {s p S a + CpCa)] 2H 3 (-1 + AH) 



+ Li 2 ((1- Vl^Ah)/2 



c 2 a A h + spc a (sf3C a - s a c )Y h 



Vl-Ah(A-h) 2 h 2 
[-2s a c aS(3Cf3 (l + s%) + sj(l + c 2 ) + S J(1 - 2s 2 )] 2h 3 (l - 40 
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+ Li 2 ((1- Vl-4A)/2) 



2c r - 



VI -AA(1 - AA)A 2 (A - h) 2 (A - H) 2 
s a c a spP 4 - [2s a c a s 3 + s}c p (l - 24)] P 3 + c p (A - #) 2 P 5 + c p s 2 a P' A 

+ Li2 I 1 " h) {A- H) 2 H 2 s l[ A2 (-3 + SH)-AA(-l + H)H-H 2 ] 
+ [4(1 + 4) - 4(1 - 24) + 2s pC ps a c a (l + 4)] 2# 3 
+s /3 s Q (s /3 s a + C/3 c a )2iJ(A(l - AH) - H) 



- Li, 1 - 



1\ (/i-l) 5 



hj (A- h) 2 h 2 
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- Li, 1 - 



C/3 



a; A 2 (A-/i) 2 (A-//) ; 



s a c a s p 2(A - l)A(h - H) 



2A A - 7 A 3 



+A 2 (1 + 5h + 5H- 2hH) -A(h + H + 3hH) + hH 



+ [2s Q c Q 4 + 4 C /3( 1 - 2s D] 2 ( A - ±)A 2 {h -H) A 3 + A 2 {1 -2h- 2H) 



+3AhH - hH 



+ c p {A - H) 2 A\-3 + 6h) - AA 3 (l + 3h + h 2 ) 
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+ log 2 (h) 
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+ log 2 (H) 
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+ log 2 (A) 



C/3 



2(1 - AA) 2 A 2 (A - h) 2 (A - Hf 



a c a Vl-AAP 4 
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+4 [-2sas a c a - Cj9 (l - 24)] Vl-4AP 3 + ^(A - if) 2 2(1 - 4A) 2 (A - hf 



+ log(/i) 



1 



2/i 2 (A - /i) 2 



[-2 W (1 + s})s a c a + s}c 2 a + 4(1 + s 2 (l - 2s 2 ))] Ah 4 



+s f3 c a (c p s a - s p c a )Ah 2 (A(l + Ah) -h) + c 2 a A' h 



+ log(H) 



2H 2 (A - H) 2 



s 2 A + 



si + sps a (sps a + c p c a ) - si(l - 2s 2 ) 



+s pCp s a c a {l + 2s 2 ) 4ff 4 + sps a {sps a + cpc a )AH 2 (-A(l + 4//) + if) 



+ log(A) 



(A -if) 2 A 4 (-76 + 48/i) 



2A 2 (-1 + 4A)(A - /i) 2 (A - H) 2 

+A 3 (32 + 156/i - 32h 2 ) + A 2 (-3 - 68/1 - 8Ah 2 ) + 6A/i(l + 6h) - 3h 2 
+s 2 4A 2 (-l + AA)(h - H) [3A 3 + A 2 (l -2h- 2H) — A(h + if - hH) + hH] 



+s 



4A 6 (-19 + 8h + AH) - AA 5 (-8 - AOh - 37H + 2AhH + 8H 2 ) 

+A 4 (-3 - 6AH - 68h - 68H 2 - 92h 2 - 312hH + 16h 2 H + 96hH 2 ) 
+A 3 (QH + Qh + 32H 2 + 3Qh 2 + 136/iff + 180h 2 H + UAhH 2 - 32h 2 H 2 ) 
+A 2 (-3h 2 - 3H 2 - \2hE - 72h 2 H - Q8hH 2 - 8Ah 2 H 2 ) + QAhH(h + H + QhH) 



-3h 2 H 2 



+ sls 2 a {-AA 2 ){-l + AA)(h -H) A 3 + A 2 {1 + 2h + 2H) 



A(h + H + hhH) + hH + sj(l - 2s 2 )4A 3 (-l + AA)(h - H) A 2 - 2A(h + H) 



+3hH 



+ spc p s a c a {-AA 2 ){-l + AA)(h - H) [2A 3 + A 2 - A(h + H + 2hH) + hH] 



+s 3 3 cps a c a (-8A 3 )(-l + AA)(h - H) [A 2 - 2A(h + H) + 3hH] 



+ lo 



(1 - Vl-Ah)/2 



- c 2 a A h + 



-2s, 



Vl-Ah(A-h) 2 h 2 
+S/3S a (sf3S a + cpc a ) - sj(l - 2s 2 ) + spc p s a c a {l + 2s 2 p ) 2h 3 {-l + Ah) 

+ [sf3S a {S(3S a + C(3C a ) - s|] Y h 



- log 



(1- Vl-4ff )/2 



VI -AH(A-H) 2 H 2 



slA H + 



4(i + 4) 
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-4(1 - 2s 2 a ) + 2spCps a c a (l + si) 2H\-l + AH) + s p s a (s p s a + c p c a )Y H 



- log (1 - Vl -4A)/2 



C /3 



1 - VI -4/A 4V1 -4/i(l -4/i) 



J (1 - 4A)Vl-4AA 2 (A - /i) 2 (A - H) 2 
+4 [-2 S/3 s a c a - C/3 (l - 2sl)} P 3 + c^(A - H) 2 P 5 + c^Pj 

+ log 
+ log 
+ log 



SpS a C a P A 



1 + y/l-Ah J 2h 2 
l-y/l-AH \ s 2 a Vl-AH(l-AH) 



1 + y/l-AH J 2H 2 
1 - y/l-4A \ c 2 ^\-AA 



1 + y/l-AA J 2A 2 
■c 2 (A-l) 2 



+ log(| - l + l/AI) 
+ 



A 2 



x 



6(1 - AA) 2 A 2 h 2 H 2 (A - h) 2 (A - H) 2 (-l + Ah)(-1 + AH) 
6(1 - AA) 2 AhH(A -h)(A- H)(-l + 4/i)(-l + 4//) x 

if(^4 - if) [A 2 (-4 + 25h) - Ah(7 + 23h) + llh 2 ] + 4(-2A)(/i - H) 

[2A 2 + A(-2H -2h + hH) + 2hH] + s\hH \a 2 {\\ -Ah + 2H) 
+A(-11H -llh + 2hH) + + s}s 2 a Hh(QA 2 )(h - H) 

+st(l - 2s 2 a )(2A 2 hH(h - H)) + s p ds a c a AA 2 hH{h - H) 



+ 7T 2 ((1 - AA) 2 (1 - Ah)(l - AH)(A - hf(A - H) 2 h 2 H 2 ^2A 2 {h - H) [(-1 + 2s 2 )4 



+2s a c a sp(l + s 2 ) C/3 J + 2 + 4 [-2 + A 2 {3 + 2h{-2 + s 2 a ) - 2H{\ + s 2 ) 
+ (1 - AA) 2 A 2 h 2 {A - h) 2 {-l + 4/i)Vl-4/f{4A H + s p s a (s p s a + c p c a )Y H 
+ [4(1 + si) - 4(1 - 2s 2 ) + 2 S/3C/3Sa c a (l + 4)] 2P 3 (-1 + 4//)} 
-(1 - 4A) 2 A 2 P 2 (A - P) 2 (-l + AH)VT Zr Ah{ - c 2 a A h + s p c a (-s p c a + c p s a )Y h 
+ [44 - 24 - 4(1 - 24) + 2spcps a c a (l + 4)] 2/i 3 (-l + Ah)} 



Vl-AAh 2 H 2 {l - Ah)(l - AH)j {(A- H) 2 P 5 + s\P[ + s 



— AA 7 
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+A 6 (-54 - 152/i + 24H) + 2A 5 (17 + lUh + 48H + 46h 2 + 6H 2 + 136/iif) 
-A 4 (5 + 96/i + 66H + U2h 2 + 58H 2 + 432hH + 168h 2 H + 184hH 2 ) 
+2A 3 (Qh + 5H + 29h 2 + 17H 2 + 94hH + 136h 2 H + 118hH 2 + h4h 2 H 2 ) 
~A 2 (7h 2 + 5H 2 + 24hH + lUh 2 H + 9QhH 2 + U6h 2 H 2 ) 



+2AhH(6H + 7h + 29hH) - 7li 2 H 2 + sls 2 a 2A(-l + 4A)(h - H) 



22A 4 



-A 3 (ll + lOh + 10H) + A 2 (l + 8h + 8H- 2hH) - A(h + H + 5hH) + hH 
+s 3 p(sp - 2sf3S 2 a - 2cf3S a c a )P 3 + s /3 c^s Q c Q P 4 |^ J. 



A. 2 The result for the 0(a%), G(a t ab) and 0{a\) corrections 

In the following we list the full result for the 0(a 2 ), 0(a t ab) and 0(a b ) corrections. As 



explained in Sect. pT4| , it has been obtained by using the Higgs sector relations eqs. ( p4[) and 
fl25|). We give this result for arbitrary space-time dimension D, using the shorthands 



i + jD + kD 2 + lD 3 



A 



i + jD + kD* 



(37) 



The result is expressed in terms of the one-loop scalar integrals Ao(m) and Bo(p 2 , m±, m-i) 
as defined in Ref. || and the scalar two-loop vacuum integral Ti 34 (m^, m 2 ,, ml) as defined in 



Ref. [221. 



SUSY,a? 



3 G% 



Higgs (4tt)W 
A (m t ) 1 8 (D - 2) m 2 B (m 2 b , 0, m t ) 
16 (D - 2) m 2 {M\ - 2 ml) m 2 s 2 B (m 2 b , M A , m b ) 8 (D — 2) m 



2^ 



Cp [m b - m t 



2\ 2 «2 



+ 



c 2 (m 2 - m 2 ) 

4 m 4 t + mt 4 + m\ m 2 (l + 2 s\ - 2 sj) - M 2 A ( C J m 2 t + m 2 b s^) ) B (m 2 b , M A , m t ) 

8(D-2) ml (4 m 2 b - w|) m 2 B (m 2 b , m b , m h ) 
(m 2 - m 2 ) 2 

(D - 2) Re(B (m 2 , 0, m b )) {{D - 2) D m\ - 2 (D - 4) D m 2 m 2 + m\ A 8 ,_ 6jl ) 



mi 



(D-2) 



(3* 



cjm^m^j + m^m? (l + 2sJ-24) - M 2 (c%m 2 t +m 2 b s% 



c 2 p (m 2 b m t - m 3 ) 2 s 2 L 
Re( J B (m 2 , M Aj m 6 )) ({D - 2) D m 4 b - 2 (D - 4) D m 2 b m 2 t + mt A 8 ,_ 6)1 ) 
2c 2 (D-2) (M\-2m 2 ) B (m 2 ,M A ,m t ) 



+ 



(m 2 b - m 2 ) 2 s 2 



■((D - 2) Dmt - 2(D - 4) Dm 2 b m 2 
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((D - 2) D mi - 2 (D - 4) D m 2 m\ + m\ A 8i _ 6>1 ) 



(-48m 4 (m 2 -4m 2 ) 



+ 



+ - 4 ° (m " ) 4 ° (m,) ( m „-2 mt ) (m» + 2m,) (-m2 + m ?) 
+ m 2 (192 m A t +Am 2 h m 2 t A_ 40 , 24 -13,2 - m£ A_ 24 ,i6 -8,1) 

+ m 2 (-384 ml- Am 2 h m t A_ 20 ,i6 -9,2 + m£ A 8 , 8 ,-4,1) ) 

2A (m f ) 2 

c 2 (-3 + -D) (Ma - 2 mt) (M4 + 2 m t ) (-m h + 2 m t ) (m h + 2 m t ) (-m 2 + m 2 ) s 2 

-2m 2 \{D-2)m 2 b \ -Am 2 (s 2 (-3D 3 + A_58.17.10) + 4 A 22,-i7,3 + A 12i8) _ 13 , 3 ) 



+ m 2 (4 (-4D 3 + A_ 52i6il6 ) + A 12i8i _ 13 , 3 + 4 A i6,-6,-3,i) 
- cl m 2 (m 2 (3 -D 4 + 4 (5 D 3 - D A + A 128) _ 10 8,i 6 ) + A_ 72i60 , 30 ,_ 23 ) 
+ 4 m 2 (-3 -D 4 + 4 (3 D 3 + A_ 140 , 136 ,_ 39 ) + A 72i _ 60 ,- 30 , 23 ) ) j 



+ M 2 



|m 2 ( - c 2 m 2 (D 4 + s 2 (2 D 3 + A 52 _ 32 _i) + A_ 24)8)24 ,_n) 



+ (D-2) m\ ((-10 + 3 D) + D A 6 ,_ 5)1 - s 2 A 2i5) _ 5) i) ) 

+ 2 m 2 ((£> - 2) m 2 (-2 D A 6 _ 5) i + s} A 10 ,-i,_ 4 ,i + sj A 14i5 ,- 6 ,i) 

-m\ (s| (_D 4 + A-164,108,27,-18) "2 (Z) 4 + A_24,8,24,-ll) + 4 + A 116i _ 92 , 21 ,_ 4 ) ) 



+ 



2A (m 6 ) : 



c 2 (-3 + D) (M A - 2 m b ) (M A + 2 m 6 ) (2 m fe - m fc ) (2 m b + m fc ) (m 2 - m 2 ) s 2 



M\ 



2m\ 4 (-3D 4 + s 2 (-4 J D 3 + J D 4 + A 116 ,-92,2i) +A_ 68i76i _ 69i26 ) 

- (D - 2) m 2 m 2 (4 (-3 + D) + (-10 + 3 D) 4 + 4 A 22,-i,-5,i) 
+ m 2 (m 2 4 (.D 4 + 4 (-2 D 3 + A_ 52 , 32il ) + A 28i _ 24 , 23i _ 9 ) 

- 2 (D - 2) m 2 (-8 (-3 + D) + sj A 14 , 5 ,- 6 ,i - 4 Aas.g.-ie.s) 



+ 2mf 



(D - 2) m 2 h m\ (8 (-3 + D) + 2 s 2 A 10)3 ,_ 5 ,i + 4 A 16 ,_ 6 ,_ 3il ) 
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+ 4m 4 b s 2 ( 3D4 + 4 ( 3 ^ 3 + A -140,136,-39) + ^68, -76,69,-26] 

+ m 2 ( - 4 (£> - 2) m\ (8 (-3 + L>) + s% A 22 ,_ 17 , 3 + sj A 14 , 17 ,_ 16 , 3 ) 

+ m 2 h sl (-2 (_D 4 + A 2 8,-24,23,-9) + s| (~-D 4 + Ai28,-108,16,5)) ) 



+ A (m 6 ) ^ (£> - 2) £> 



2-D 4 ((-l + J) m 2 + m 2 ) 



4m^m 2 - 



(m 2 b -m 2 t f 
(M\ -m\-m\) (c 4 m 2 + m 2 S J) 



r 2 2 



B (m b ,M A ,m t ) 



, 2 8(D-2)m 2 b / 4 4 , 4 4 



+ 8 (D - 2) m 2 . Re(B (mj, 0, m 6 )) + 



C/3 m t + m fe ^ 



4 ( m l - m ?) 2 4 

+ m 2 m\ (l + 24-2sJ)- M 2 (cj m 2 + m 2 4) ) Re(5 (m 2 , M A , m b )) 
16 4 (J - 2) m 2 m 2 (Mj - 2 m 2 ) E (m 2 , Ma, m t ) 
(m 2 - m 2 ) 2 4 
_ 8(D -2)m 2 b m 2 (m 2 h -Am 2 t ) B (m 2 ,m h ,m t ) 
(m 2 - m 2 t f 

(D - 2) B (m 2 , 0, mt) (-2 (L> - 4) D m 2 m 2 + (D - 2) D m\ + m 4 A 8 ,_ 6)1 ) 



+ 



2 (-D — 2) (M 2 - 2 m 2 ) 4 ff (m 2 , M A , m 6 ) / 2 2 

— — 2 [-2(D-4)Dm b m t 

4 (m 2 b -m 2 ) V 

+ (D — 2) Dmf + ml A 8 ,_ 6il ) - (D = 2) 

' (mg - m 2 ) 

( - 2 (0 - 4) D ml ml + (D - 2) D mf + m t A 8 ,_ w ) - m , (4 ^^ff w _^ 
(l92 + 4 m 4 (-96 m 2 + m\ A_ 40i24i _ 13i2 ) + (192 m\ - m 4 A_ 24 ,i 6 ,_ 8 ,i 

- 4 m 2 m 2 A_ 20) i6 _ 9)2 ) + m 2 h m\ (-48 m\ + m 2 A 8i8 _ 4ji ; 

A (m t ) 



+ 



44 m 6 m t (-^1 ~ ( m < ~ m b) 2 ) 2 (^A - ( m < + m fe) 2 ) 2 ( m 6 - m t , 



{m 2 b - m 2 ) 3 ^4 ( D ~ 2 ) 2 Dm S t +(D- 2) 2 Dm\ s} + m 4 m 4 ( - (10 - 3 D) 2 D 
+ 96 (-D — 2) 4 - 96 (£> - 2) sj) + m 2 ((£> - 2) 2 D + 48 (£> - 2) S J - 2 s 2 A_ 48i60i _ 32 , 5 ) 
+ m 2 m\ (48 (D - 2) s% + L> A_ 68i60 ,_ 9 + 24 A 48)12) _ 32 , 5 ) 
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+ M\ ^-c 2 p (D-2) 2 Dm 6 t + (D-2) 2 Dm 6 b s 2 p + m 2 b ( - 32 (-1 + £>) + s£ A_ 

- 4 sj A_ 32i 4o -18,3) + ml m\ ((£> - 2) 2 D + 4 sj A_ 32i40 _ 18j3 + sj A 96 _i 32 , 76 -13) 

-AM\(^-cl{D-2) 2 Dml + {D-2) 2 Dmlsl-Amtmt (-1 + 2^) A 8 ,_ 10 ,i 
+ m 2 (-4 A_4 A i + 4 A_i 28 ,i48 -52,9 - 2 sj A_ 56i68 _ 30 , 5 ) 

+ m b m t (iP — 2) 2 D + 2Sp A_5 6!68 _ 30 ,5 + 4 A 96 -124,68 -ll) ^ 

-4Mi(m2-m?) 2 ^-cJ( J D-2) 2 J Dm« + ( J D-2) 2 J Dm«4 

+ 4mtm 4 t (-1 + 2 s 2 ) A_ 4 ,i 8 ,-8,i + ™b rn 2 ((D - 2) 2 D + 4 sj A_ 24 , 2 o,-6,i + 4 A 9 6,-ii6,60,- 9 ) 
+ m 2 (4 L> A 8 _ 8) i - 4 sj A_ 24i20 _ 6 ,i - 4 A 96 _ 44 _ 12> i) 

+ 2M^-3c 2 (D-2) 2 Dmf + 3(D-2) 2 Dmi°s 2 -2m^m^ (sj (17D 3 + A_ 72j204i _ 118 ) 

- 4s 2 (5D 3 + A_ 32j78i _ 36 ) + A_ 40 , 84 ,-26,3) +2mtm 6 t (l6 -2AD + 4 A_ 72i204 ,-ii8,i7 

- 2 S 2 A„ 8i48i _ 46 j) + m 2 m t 8 (s 2 (13 -D 3 + A_ 320 , 29 2,-68) + 4 (~ 22 ^ + A 304,-328,13 2 ) 

+ Ai 6i24i _5 2j6 j + ml m 2 (3 (D — 2) 2 D + 4 A_ 304 , 328 , -132,22 + 4 A 288j _ 364i i9 6i _ 3 i) 



+ A (M A ) < - A (m t ) * ^ 

[ c 2 (Ma - 2 m t ) (M A + m 6 - m t ) 2 mf (Ma - m fe + m t ) 2 

1 

^-M A + m 6 + m t ) 2 (Ma + m b + m t f (M A + 2 m t ) (-m 2 b + m\) s\ 
M\° (^-(D-2) 2 D mt sj + m 2 b m 2 (-3 (D - 2) 2 D + 6 (D - 2) 2 £> s£ - 2 sj A_ 24i20i _ 2>1 ) 

+ 4 m A t A 16il2 ,_ 24)3 j +2M\( [ 2{D-2) 2 Dmtst ) + m 4 b m 2 t (D — 2) 2 D — 10 (D — 2) 2 D si 
+ 4^72,68,-14,7) + m 2 < (4 (6D 3 + A_ 192)136il6 ) -44 (5 D 3 + A_ 4j8 _ 5 ) + A_ 24;12 _ 2;9 ) 

- cj m\ A 88 , 20 ,-94,i5 j + 4 m 2 (m 2 b - m 2 ) 2 ^(D - 2) 2 D m\ sj + m£ m 2 (2 (D - 2) 2 (1 + 2 D) 

- 4 4 A 8 _ 2 _ 3) i - 4 ^8,52,-45,5) - 2 4 ™t A 20 -4,-13,3 + - 4 4 A_ 52) n 2 _ 45)3 
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+ sj A_ 24 ,100 -54,1 - 2 A 44i 20 -39,7) + ™fc ™t (sj (9 D 3 + A 72 -60 -14) +4sJ A-64,114 -36,1 

+ A 12 o,-36o,i66,-i 6 )j -2 Ml ^3(D-2) 2 Z}m*4 + m 2 m t 6 ( S J (17 D 3 + A_ 648j588 , 2 ) 

+ (-44 -D 3 + A 256 ,_ 288i8 ) + A_ 192i36 , 92 , 9 ) + ml m 2 (9 (D - 2) 2 D - 12 (D - 2) 2 L> S J 
+ 4 A_ 88i68i _ 10i h j + m\ (s^ (17 -D 3 + A_ 248jl08i78 ) — 4 4 

(10D 3 + A 20 ,_ 16i _ 5 ) 

+ A_ 24i32 ,- 46 ,3o) - 2 cj m* A 148i4 ,-i3i,24 j + 2 M^2 (7J - 2) 2 /J m* S J + m\ m\ (4 s 2 (l9 D 3 

+ Ai3 2i 6 4j „i 4 3 ) - 2 sj (37 D 3 + A 168i236i _ 298 ) + A_ 328i84 ,ii4,3) + m\ m 2 (6 (-D - 2) 2 /J 

- 6 (/J - 2) 2 D + 4 A_ 56i20i6i7 ) + m 2 m t 8 (2 4 (49 D 3 + A_ 304 , 612 ,_ 248 ) 

- 4 4 (38 D 3 + A 44il52i „ 139 ) + A_ 8i _ 164i82i25 ) + ml m\ (4 4 (9 /J 3 + A_ 24i52i _ 14 ) 

- 4 4 (14 D 3 + A 52 ,l 6 ,-4l) + A 8i20 ,-66,35) - 4 m <° A 440,- 

- M\ \[D - 2) 2 D m 12 4 + 2 m£ m° (4 4 (7 /J 3 + A_ 104 , 2 32,-io 2 ) + 2 4 (59 -D 3 + A 432 , 4 ,_ 276 ) 

+ A_368,-276,408,-65) + mj° m 2 (3 (D — 2) 2 D — 2 (D — 2) 2 D + 8 sj A_ 4)2 ,_i,i) 
+ mgm* (4 (-7D 3 + A 32 ,- 268 ,244) - 24 (47 D 3 + A192.156.-212) + A 64i44i _ 124i47 ) 

+ m£m* (4 (-221 D 3 + A 1088 , -2900, 1380 ) + 4 4 A_i6, -196, 140,1 + Ai60, 1272, -1040,142) 

- 4 "lj 2 A 608 ,-92,-428,93 + ™ft ^sj (-330 L> 3 + A_ 2496>1272i968 ) 

+ 84 (32 -D 3 + A 44il50) _i53) + A 1120) -i 892i444i 3i 



-36,-342,69 



+ A (m b ) 



c 2 m 2 (M\ -Am 2 ) (M A + m b - m t y (M A - m b + m t ) 1 (-M A + m b + 



1 



M) cj (£> - 2) 2 Dmf-mt 4 A 16i12i _ 24 , 3 



(M^ + m b + mt) 2 (ml - m 2 t ) s 2 



'13 

+ m 2 b m\ (A_ 48) 28,8,-i + 24 A_ 24 , 20 ,-2,i + 2 4 A 48 ,_ 2 8,-8,i) j +A(m 3 b -m b m 2 t ) 2 

4 (- ~~ 2 ) 2 Dm 8 t +ml ml (-2 4 (5 L> 3 + A_ 80 .i 24 _ 36 ) + A_ 96i388 _ 204 ,25 - 4 A - 2 4,ioo,-54,i) 

+ 2 mg sj A 20 ,-4,-13,3 + ™6 ™t (4 ( 5 ^ 3 + A 8 ,52,-4 6 ) - 2 4 (7 D 3 + A 24 , 48i _ 52 ) + A 32 ,36,-44, 5 ) 

+ m\m\ (4 (-9D 3 + A_ 72i60 ,i 4 ) +24 (ll D 3 + A_ 56)168 _ 84 ) + A 64 _ 36 _i 2i3 ) 
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+ 2 M\ ^ - 2 c% (D - 2) 2 D m\ - m\ m\ [s% (6 D 3 + A_ 192il36 , 16 ) + 4 4 (2 D 3 + A 92 ,_ 60 ,_ 13 ) 
+ A_2oo,ii6,34 -5) + m 2 b m\ (sj (-7 L> 3 + A 72 -68,14) + 4 4 A_ 36i24 ,3,i + A 72 _ 52 _ 2 _ 3 ) 

2 M\ ^ — 3cp(D — 2) 2 Dm^— mf m\ [s% [\1 D 3 + A_ 648 , 588 , 2 ) 
+ 2 (5 L> 3 + A 520i _444,- 6 ) + A 

-584,336,102,-18 

)-mimt ( S J (l7D 3 + A 

-248,108,78) 

+ 2sJ (3 D 3 + A 288i _ 140 ,-88) + A_ 352i2 04,52,7) + m 2 b m t 6 ( - 2 (4 £> 3 + A_ 4 4,28,i) 



+ ^ si A, 



°f3 ^88,20,-94,15 



+ 24 A_ 88i44 ,i4,5 + 4 ^88,-68,10,-11 J + 2 sj Ai 48 ,4,-i3i,24 J + 2 - 2 cj (£> - 2y D m 

- m 2 b ml (sj (7 -D 3 + A_ 56i2 o, 6 ) - 4 4 (2 D 3 + A_ 28i4 , 9 ) + A_ 5 6,2o,6, 7 ) 

+ m^m^ (-44 (18D 3 + A 36il72i _ 155 ) +2 S J (37 D 3 + A 168j236i _ 298 ) + A 136i132i _ 138i _ 5 ) 

+ mi m\ (4 (4 D 3 + A_ 100 ,88,i3) - 4 4 (9 D 3 + A -24,52,-u) + A 296i _i 64i _4 2i _i5) 

-|- mjf 1 4 A440 5 — 3 6,— 342,69 ^Ife 

(4 4 (ll D 3 + A_ 348 ,46o,-io9) - 2 4 (49 D 3 + A 

-304,612,-248) 

+ A 792 ,_452,-i42,2 9 ) J+Ma ^4 ( D - 2 ) 2 D m * 2 + m &° ™* (8 4 (32 -D 3 + A44.150.-153) 

+ S /3 ( — 182 D 3 + Ai 792i _ 3672i i4 8 o) + A_i 02 4,580,188,-43 j + s /3 A-608,92,428,-93 

+ m 2 b m]° (-2 4 (7 D 3 + A_ 32jl2i _ 4 ) + A„ 32i20 ,_i 2 , 9 + 84 A_ 4)2> -i ) i) 

+ 2 ml ml (44 (7D 3 + A_ 104 ,232,-io2) -24 (87D 3 + A 16i932i _ 684 ) + A 80 , 6 60,-552,8i) 

+ m 4 b m 8 t ( 

A 32) _ 268)2 44 5 _ 7 + 4 

4 A 8 o, 2 i 2i - 228i27 — 2 Ai44 i268i _ 272i27 ) 

+ (2 4 (219 -D 3 + A_ 105 6,3292,-166o) + 4 ( " 221 ^ + A 1088,-2900,1380^ 

+ An 84i _ 2 4i 2i9 oo,-75 j 

+ T lu (ml, ml Ml) =— J - Mi (ci (D - 2) 

V b b A c 2 (M 2 A -Am 2 b ) (rn 2 b - m 2 ) 2 si I H^ 1 ; 



tJ S /3 



- 4 (-1 + D) ml 4 + m 2 b m 2 t s} (4 - 2 £> + (-6 + 5 D) 4) ) + M 2 A m 2 b [{D - 2) 
(6 - 8 4 + D 4) - 2 ml m\ s% (6 (£> - 2) + 4 ^22,-21,2) + ^ sj A 28i _ 32>3 ) 

+ m t i^t 4 A -36,48,-n,i - 2 m t 2 4 (4 - 4 D + A_i 6>30 -9,1) 
+ (-8(Z>-2)8 (-3 + D) s 2 + sj A 4l i 2l -7 1 i) j J 



4 
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Ami 



T 134 (ml ml m 2 b ) { m\ ((2 + D) m\ + (D - A) 2 m 2 t ) 



{Ami -m 2 h ) (™%-m$y 
+ m t ( m l A -28 -6,1 -2mj A 26 _5 -4,1) +2 ml A 34)15 _ 8>1 
+ 4T 134 K,m^0) f _ 12m 4 + m 2 (_i2^ + m 2 A _ 22)i9 _ 7)i) +m 2 (24m 2 + m 2 A_ 22 , 19 ,_ 7il ) 
+ 4T 134 K,m2,M|) 



c| (Ma + m fc - mt) 2 (M4 -m b + m t f {-M A + m b + m t f (M A + m b + m t ) 2 (ml - m 2 t ) 2 s 2 p 
l^2M A ° ^cj (-2 + 3D) m t 4 + (-2 + 3L>) m\s\ + m 2 h m 2 t (2-3D + 2 (-6 + 5D) s\ 

- 2 (-6 + 5 D) 4) j - 2 M\ (^4 ml A_ 2 , 7 , 2 + ml 4 A_ 2 , 7 , 2 - m 2 m\ [s% A_ 34 , 2 3,2 - 8 s 2 A_ 2 -1,1 
+ 4 A 2 _ 3i i) + m\ m 2 t ( (-52 + 54 D - 4 -D 2 ) 4 + (34 - 23 D - 2 D 2 ) 4 + A 10 ,_i 9 , 2 ) j 

+ (m 2 - m 2 ) 4 ^cj D ml A 6 _ 7>1 + £> m£ sj A 6 _ 7>1 - m 2 . (d sj A 6 _ 7>1 + 4 4 A_ 4i22i _ 9il 
+ A 8 _io4,50 ,-e) + ml m\ (-D sj A 6 _ 7jl + 2 4 A_ 8)50 _ 25 , 3 + A 8 ,io -7,1) j 

- M\ {m 2 b - m 2 ) 3 ^ - cj m° A_ 8 , 22i _ 25 , 3 + m° sj A_ 8 , 22i _ 25 , 3 + m 2 (4 (-2 + 5 D) 

-AD si A 56 ,- 2 5,3 + 4 A-8,154,-75,9) + ™ 4 m 2 (2 4 (3 -D 3 + A_ 8i42i _ 25 ) + sj (-9 -D 3 + A 8i _ 154 , 75 ) 
+ A 16)50 ,-25,3)j +M\ {ml-m 2 ) ^m 2 m^2sj (3 D 3 + A 12;32 _ 25 ) + s| (-18 L> 3 + A 24 _ 328)158 ) 

+ A_ 44il68 _ 75)9 ^ - cj m t 8 A_20,24,-33,3 + m b 4 A-20,24 -33,3 ~ ™6 ™t (2 sj (3 -D 3 + A_ 36>100 _ 29 ) 

+ 24 (3,D 3 + A 12i32i _ 25 ) + A 4i _ 96i33i _ 3 ) + m 4 b m 4 t (-1 + 2 4) A 36>112i _ 75i9 j 

-Ml ^4 (-12-19D 2 + D 3 ) m 8 t + (-12-19 D 2 + D 3 ) m\s%-m\m\ (a 4j32i _ 19)1 

- 2 4 A 36i64i _ 47i5 + 24 A 36i64i _ 47j5 j + m 2 (4 4 A_ 20i32i -7 5 i — 2 s| A_ 4i _ 8j5i i + A 44 -56,19,-1) 

+ ml m\ (A_28,88,-19,l +44 A_ 20 ,32 -7,1 + 4 Al52 -272,66 ,-e) 
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4 T 134 (m 2 ,m 2 h ,m 2 b ) 



mi 



+ 



12 ml + ((2 + D) m 2 h - 12 m 2 ) (m 2 - m 2 ) — m 2 ((14 + D) m\ + 2Arr 
8 



8 

c 2 (M A -2mt) (M A + 2m t ) (m 2 - m 2 ) 2 s 2 



T 134 (m t ,m t ,M A 



4 cj (-! + £>) 



- (D-2)mis 4 p + c}m 2 b m 2 t (2 - 3 D + (-6 + 5 D) s 2 ) 
+ M 2 A m 2 ^2c 2 p m 2 b m 2 ((22 - 21 D + 2 D 2 ) s 2 + A_10.15.-2) + ™6 {{D - 2) 2 + (D - 2) D s 

-2sp A 8 _ 6jl ) + cj A 28 ,-32,3^ + m< ^cj m t A_ 36;48 -11,1 + m\{^ - 2 s 2 p A_ 8) i 6 _ 7jl 

+ A_ 4 .i2 -7,1 + sj A4.12 -7,1) + 2 m 2 m 2 (_J A_i 6)30 _ 9 ,i + A i2 -26,9 -1) j J 



+ 



4 m 2 



(m/j - 2 m t ) (m ft + 2 m t ) (-m 2 b + m 2 ) 
+ m 2 ((D - 4) 2 m 2 - 2 m 2 A 26i _ 5i _ 4il ) + 2 A 34l i5 1 -8 1 i 



T lu (m t , m t , m h ) I (2 + D) m h + m h m t A_ 28 -6,1 
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